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Abstract. The recent experimental support for the presence of the Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) phase in the CeCoIns directed the attention towards the 
mechanisms responsible for this type of superconductivity. We investigate the FFLO 
state in a model where on-site/inter-site pairing coexists with repulsive pair hopping 
interaction. The latter interaction is interesting in that it leads to pairing with nonzero 
momentum of the Cooper pairs even in the absence of the external magnetic field (the 
so-called ^-pairing) . It turns out that depending on the strength of the pair hopping 
interaction the magnetic field can induce one of two types of the FFLO phase with 
different spatial modulations of the order parameter. It is argued that the properties 
of the FFLO phase may give information about the magnitude of the pair hopping 
interaction. We also show that ry-pairing and d-wave superconductivity may coexist 
in the FFLO state. It holds true also for superconductors which in the absence of 
magnetic field are of pure d-wave type. 
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1. Introduction 

An unconventional superconducting state with a non-zero total momentum of the 
Cooper pairs was predicted by Fulde and Ferrell [1] as well as by Larkin and Ovchinnikov 
[2] in the middle of 60's. Under particular conditions, this phase should occur at 
low temperatures and in strong magnetic fields. Due to sever requirements for the 
formation of the FFLO state, this type of superconductivity has experimentally been 
observed only recently. In the FFLO state the superconducting order parameter (OP) 
oscillates in the real space. This property, to some extend, resembles the unconventional 
superconductivity in strongly correlated systems [3], where the OP changes sign in the 
momentum space. The FFLO state has recently been analyzed in the context of heavy 
fermion systems [l[5l[6l[71[8l[9j[IOl[IIl[I^ organic superconductors 

PH ED], ultracold atoms [2H [22] and a dense nuclear matter [23, E5 EE] . Although, there 
is no direct experimental evidence for the spatial variation of the OP, suggestions for 
the future experiments have been developed in Refs. [9j ES, EZ] • 

The orbital (diamagnetic) pair breaking is a crucial mechanism that limits 
realization of the FFLO state. In the vast majority of superconducting materials it 
is a dominating pair breaking mechanism that destroys superconductivity for magnetic 
fields much weaker than the Clogston-Chandrasekhar limit (H cc ) [28| [29] . It holds 
true also for models appropriate to describe the short coherence length superconductors 
[30| I3TI [32] . The significance of the diamagnetic pair breaking is usually described 
in terms of the Maki parameter a = \/2H°2 b / ' H cc ', where i?°J 6 is the upper critical 
field calculated without the Zeeman splitting. There exist two general possibilities to 
reduce the destructive role of the orbital pair breaking. In the layered superconductors, 
formation of Landau orbits should be suppressed for magnetic fields applied parallel 
to the layers. This may explain possible observations of the FFLO state in some 
organic superconductors [T9l [20] . The role of the orbital pair breaking should also 
be limited in systems with narrow energy bands, like heavy fermion systems. The 
experimental evidence for the FFLO superconductivity in these systems seems to the 
strongest [l[3[g[3[8l[9l[T0l[IIJ[T2l 

In the context of recent investigations of the FFLO state, it is important to search 
for other mechanisms that stabilize superconductivity against the orbital pair breaking. 
Recently it has been found that superconductivity originating from repulsive pair 
hopping interaction is unique in that it is robust against this pair breaking mechanism 
[33J. The origin of this interaction may vary in different systems and therefore we 
do not specify a particular superconductor for which the following qualitative analysis 
can directly be applied. The repulsive pair hopping interaction can be derived from a 
general microscopic tight-binding Hamiltonian [31], but in this case the magnitude of 
the interaction is very small. However, since this interaction leads to superconductivity 
that is almost unaffected by the orbital effects, it may become more important close 
to the upper critical field, i.e., in the regime where the FFLO phase is expected to 
occur. We may also consider other sources of the pair hopping which give rise to much 
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larger magnitudes of this interaction. For example, such term may be included in the 
effective Hamiltonian describing Fermi gas in optical lattice in the strong interaction 
regime [35j EE]. It also arises in a natural way in multiorbital models [37], though then 
the pairs hop between different orbitals on the same site. Nevertheless, we expect that 
some of our conclusions can still be valid. The role of the pair hopping interaction in 
a multiorbital models is of particular interest because of its presence in the recently 
discovered iron-pnictides [38J. Additionally, if the on-site repulsion exceeds the gap 
between the lowest and the next-lowest bands in optical lattice, then the interband 
pair hopping seems to important. Very recent Quantum Monte Carlo calculations for 
TMTSF-salt [39] also suggest a significant role of the pair hopping processes in this 
system, which, on the other hand, probably exhibits the FFLO phase at high field 
[201110]. 

In the absence of magnetic field, the pair hopping interaction is responsible for the 
77-type pairing where the total momentum of the paired electrons is Q = (tt, 7r) and the 
phase of superconducting order parameter alters from one site to the neighboring one 
[HI H21 HH HH H5] . It has been shown that flux quantization and Meissner effect appear 
in this state [HI H7J. Therefore, even in the absence of external magnetic field, the 
repulsive pair hopping interaction favors pairing with non-zero momentum of Cooper 
pairs. Although, the pair hopping may not be the dominating pairing mechanism, 
its presence may affect the FFLO phase. In the following, the superconducting 
state with zero total momentum of the Cooper pair will be referred to as the BCS 
superconductivity. In that sense, both FFLO and 77-pairing will be considered as a 
non-BCS state. 

In the present paper, we analyze the role of the repulsive pair hopping interaction for 
the FFLO state. The above mentioned systems, where the pair hopping interaction may 
play a significant role, exhibit both inter-site pairing (CeCoIn 5 , TMTSF) as well as on- 
site pairing (optical lattice, the iron-pnictides). Therefore, we consider a model where 
this interaction coexists with on-site or inter-site pairing potential. However, we do not 
refer to any particular system. As the microscopic mechanism of superconductivity in 
most of the unconventional superconductors is still under debate, we do not discuss the 
origin of the pairing potentials. 

For the system with on-site pairing, magnetic field reduces the total momentum of 
electrons forming Cooper pairs in the ^-pairing state. As a consequence, the amplitude 
of the superconducting order parameter becomes a site dependent quantity. Recent 
theoretical investigations of the FFLO phase in the attractive Hubbard model have been 
motivated mostly by the increasing interest in the ultracold Fermi gases [48j|49]. These 
approaches may also be applicable to compounds other than the strongly correlated 
heavy fermion systems [501 EE]- Contrary to this, experimental results obtained for 
the CeCoIns indicate on the anisotropic d-wave pairing. Therefore in Sec. IV. we 
study the case of inter-site attractive interaction that is responsible for the ci-wave 
superconductivity. We show for the case of inter-site attraction that <i-wave and re- 
pairing orders may coexist in the FFLO state. Such a coexistence is possible also for 
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systems which in the absence of magnetic field are in the pure <i-wave state. For on-site 
and inter-site pairings, the potential of the pair hopping interaction J is assumed to be 
positive. 



2. On— site pairing 

We start our analysis with a model with on-site pairing interaction described by the 
following tight-binding Hamiltonian: 



h= -t c I* c 3<t + u 4^4^! 

+ J 4Ai c n c n ~ J2 \. s ( a ) h + fA C L C ^ C 1 ) 

{id) i> a 

where t is the nearest neighbor hopping integral, J is the pair hopping interaction, \i is 
the chemical potential, s(f) = 1, and s(|) = —1. The Zeeman coupling is determined 
by h = gfiB7~L/2, where g is the gyromagnetic ratio, \lb is the Bohr magneton and H 
is the external magnetic field. Hamiltonian ([1]) does not include the diamagnetic pair 
breaking. We refer to Refs. [52] ESI El EHl M> EHJ ESI EDI EE] for the discussion 
of the influence of this mechanism on the FFLO state. Here, we focus on the role of 
the pair hopping interaction as well as on the properties of the ^-pairing, that is robust 
against the orbital pair breaking. Nevertheless, the role of the orbital pair breaking will 
be briefly discussed. In this section we assume the simplest form of the effective on-site 
pairing interaction {U < 0), that is responsible for the s-wave superconductivity. 

We apply the mean field approximation and assume the order parameter in a general 
form: 

M 

A(Rj) = (cjiCtf) = A m exp(iQ m • Rj). (2) 

771=1 

Then, the Hamiltonian in the momentum space takes the form 
Hmf = X! £fco- c L- c fc°- 

k<i 
M 

+ U es(Qm) X) ( A m C -fc+Q m |Cfc T + H.C.) 

m=l fc 

M 

- N U cS (Q m )\A m \ 2 (3) 



m=l 



where 



£ka = £k - fJ> - s(a)h. (4) 

For arbitrary lattice geometry the dispersion relation is given by = 
~ t/NYlj exp (iRj • k), where the prime means summation over the nearest neighbor 
sites. We have also introduced an effective pairing potential 

U eS (Q) = U-^-. (5) 



FFLO phase in the presence of pair hopping interaction 



5 



For a general form of the order parameter [Eq. §2§] diagonalization of the mean 
field Hamiltonian usually cannot be reduced to an eigenproblem of a finite Hermitian 
matrix. Therefore, we restrict further discussion to two simplest cases. The first one 
was originally proposed by Fulde and Ferrel (FF), whereas the second by Larkin and 
Ovchinnikov (LO). In the former one, it is assumed that M — 1, so the absolute value 
of A(Rj) is constant, but the phase changes from site to site. In the latter case M — 2, 
Ai = A 2 and Qi = —Q2- Then, one gets A(Rj) = 2A cos(Q ■ Rj), where we use 
Ao = Ai and Q = Q±. However, one should keep in mind that at low temperature and 
high magnetic field FFLO phases with M > 2 may be thermodynamically more stable 
[621 153] ■ 

For the FF phase, the mean field Hamiltonian can be diagonalized by means of the 
Bogoliubov transformation. Straightforward calculations lead to the following form of 
the grand canonical potential Q = — ZcTlnTr exp(— /3H) 

n = - kT J2 E ln I 1 + eM-PE hta )} 

a=± fc 

+ J2^+Ql-NU eS (Q)\A \ 2 , (6) 



where 



Ek,± — 2 [ £fc T _ e -k+Q[ 



± + £-k+Qi) + 4f/ c fr(Q) 2 |A 



o 



12 



(7) 



In the case of LO superconductivity, the Hamiltonian can not be diagonalized 
analytically. However, the pairing term links the one-particle states with momenta 
lying along a single line in the Brillouin zone. Therefore, one can solve the resulting 
eigenproblem numerically for relatively large clusters. Namely, for a L x L cluster, one 
has to diagonalize a 21 x 21 Hermitian matrix. 



2.1. Numerical results 

We start our discussion with the simplest case M = 1 (FF state), that allows one to 
estimate the boundaries of the non-BCS superconducting phases. However, the presence 
of the LO superconductivity will be discussed as well. The thermodynamically stable 
phase has been determined through minimization of the grand canonical potential with 
respect to the superconducting order parameter |A | and Q. The calculations have 
been carried out for a square lattice with \x = as well as for a triangular lattice with 
fi = 2t. These chemical potentials correspond to maxima in the density of states and, 
therefore, to the highest superconducting transition temperatures. A comparison of 
results obtained for both the cases allows one to check the role of the lattice geometry. 

First, we focus on the influence of the pair hopping interaction on the properties of 
the FFLO phase. Fig. [[] shows how the ground state of the system depends on J and 
h. In Figs. [2] (square lattice) and [3] (triangular lattice) we present the values of Q that 
minimize fl for some particular values of J and h. In the absence of the magnetic field 
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Figure 1. Phase diagram showing the stable superconducting phases for various 
J and h at T = 0. The upper panel shows results obtained for a square lattice with 
U = —2.0t, whereas the lower one shows results for a triangular lattice with U — —2.5t. 
The dashed lines are explained in the text. 

there are two stable superconducting phases for both the lattice geometries. For small J 
there exists an isotropic BCS phase, that will be referred to as s-wave superconductivity 
(see Figs. [2^ and[3k). The ^-pairing phase occurs for larger J. In the case of a bipartite 
square lattice, ^-pairing corresponds to the total momentum of Cooper pairs Q = (n, n), 
where the phase of superconducting order parameter changes from one lattice site to the 
neighboring one (see Fig. [2b). Stability of this phase obviously follows from the fact, 
that the pair hopping interaction involves sites, which belong to different sublattices. 
Then, the oscillating Aj minimizes the energy of the system, for the physically relevant 
repulsive interaction. In this context, the problem of //-pairing on a non-bipartite 
triangular lattice is interesting even in the absence of magnetic field [64J. We have 
found that also for this geometry, the repulsive pair hopping interaction may lead to 
thermodynamically stable phase with Q ^ 0. The ground state energy is minimal when 
Q represents one of the corners of the hexagonal first Brillouin zone (see Fig. [3b). It is 
easy to check that for such a value of Q, the phase of Aj takes on three different values, 
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Figure 2. Grand canonical potential minimized with respect to |Ao| for given 
Q = (QxiQy)- We have used the same parameters \x and U as in Fig. [TJ whereas 
the values of h and J are shown explicitly above the small panels. Note different 
ranges of Q x and Q y in various panels. These results have been obtained for a square 
lattice. 



namely, Aj = A , Aj = A exp(i|7r) or Aj = A exp(— depending on i. In the 
presence of sufficiently strong magnetic field there exist four different superconducting 
phases. Apart from the discussed above s-wave and 77-pairing states there are other 
two phases, that will be referred to as FFLOl and FFL02. Investigation of the total 
momenta of Cooper pairs (see panels "b" and "d" in Figs. [2] and [3]) , allows one to link 
FFLOl and FFL02 to s-wave and 77-pairing, respect what, in turn, is a hallmarkively. 
Namely, Q obtained for the FFLOl is relatively close to the origin of the Brillouin zone, 
whereas in the FFL02 phase Q remains on the edges of the zone. One can see, that 
FFL02 phase occurs for lower magnetic fields than the FFLOl. FFLOl evolves from 
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Figure 3. The same as in Fig. [5] but for a triangular lattice. 

the s-wave superconductivity when sufficiently strong magnetic field is applied. It holds 
true for both the lattice geometries. On a square lattice, FFL02 evolves from the re- 
pairing state under the same conditions. However, on a triangular lattice this phase is 
stable only for moderate values of the pair hopping interactions as well as for moderate 
magnetic fields. It is surprising that increasing of magnetic field may cause two phase 
transitions, the first one from s-wave to FFL02 is discontinuous and the second from 
FFL02 to r/-pairing is a continuous transition (see Fig. H]). We have also found that 
the transitions from the s-wave phase to the FFLOl state are discontinuous, whereas 
the transitions from the FFLOl and FFL02 phases to the normal state are continuous. 
Finally, for a square lattice the transition from the r/-pairing state to the FFL02 phase 
is discontinuous as well. The field dependence of the superconducting order parameter 
in various phases can inferred from Fig. HI 

In the absence of the pair hopping interactions, it is well known that the LO phase 
has lower ground state energy than the FF one [2]. In our case, this result directly 
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Figure 4. |t/ e fr(Q)Ao| as a function of magnetic field. The parameters as well as the 
lattice geometry are explicitly shown in the figure. 



applies to the FFLOl phase. We have found that in the case of FFL02, state with 
M = 2 has energy lower than that with M — 1. Therefore, also in this case |Aj| 
is spatially inhomogeneous. This result has been obtained on the basis of numerical 
diagonalization of 200 x 200 clusters with periodic boundary conditions. Consequently, 
there exists a simple criterion to distinguish between FFLOl and FFL02 phases. In 
both the cases A* ~ cosQ • Ri but \Q\ < 1 for FFLOl, whereas \Q\ > tt for FFL02. 
Therefore, the periods of spatial modulations of the order parameters relevant to FFLOl 
and FFL02 are very different. In the latter case it is of the order of the lattice 
constant. It is instructive to examine in more detail the spatial modulation of Aj 
in the FFL02 phase on a square lattice. As the total momentum of Cooper pairs Q 
is close to II = (vr,7r), one can introduce Q' = II — Q and note that \Q'\ <C 1. Then, 
Aj ~ cos[(II— Q')-R^ = cos(U ■ Ri) cos(Q' -Ri). The spatial profile of Aj is determined 
by two oscillating functions. Due to the first one the phase of the superconducting order 
parameter alters from one site to the neighboring one. It means that the FFL02 phase 
retains the basic properties of the ^-pairing superconductivity. The second factor is 
responsible for a slow variation of the magnitude of superconducting order parameter 
| Aj| what is a hallmark of the LO-type of superconductivity. On the basis of our analysis 
one can not exclude that the FFL02 phase with M > 2 is more stable. Therefore, the 
actual boundaries of the FFL02 state may cover a slightly wider range of magnetic 
fields, than presented in Fig. CD 

Up to this point we have analyzed a two-dimensional system, where the influence of 
the orbital pair breaking can be neglected provided the applied magnetic field is parallel 
to the plane. However, if the magnetic field has a nonzero component perpendicular 
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Figure 5. |Ao| as a function of magnetic field and temperature for a square lattice 
with J = 0.35*. 



to the plane, as well as in the case of three dimensional systems the destructive role 
of the orbital effects has to be taken into account. This problem has been analyzed, 
e.g., in Refs. §Q\ EH E2j [331 E2j E31 E2 ESI ESI E3 ESI E31 EQ1 EI] ■ On the one hand, it 
is the dominating pair-breaking mechanism for s-wave and FFLOl superconductivity 
[301[ID3[33E21E31E1E^ On the other hand, it is very ineffective 

in destroying 77-type superconductivity [33]. Therefore, we assume that FFL02 will 
also be robust against the orbital pair breaking. This property may lead to significant 
modifications of the phase diagram presented in Fig. [TJ Namely, the dashed lines show 
the boundaries of the rj and FFL02 types of superconductivity obtained under the 
assumption that s-wave and FFLOl phase are destroyed by orbital effects. Note, that 
even weak pair hopping interaction should lead to the onset of rj and FFL02 phases for 
fields sufficiently strong to destroy the conventional superconductivity. However, this 
conjectural result should be confirmed by calculations for the FFL02 phase with the 
diamagnetic pair breaking explicitly taken into account. 

Finally, we discuss the standard (kT, h) phase diagram for a square lattice with 
J = 0.351 It is the value of pair hopping interaction for which the ground state in the 
absence of magnetic field is the r^-type superconductivity (see Fig. [[]). Fig. E] shows the 
results. It is interesting that despite unconventional character of the 7/-pairing, the phase 
diagram is very similar to analogous phase diagram for BCS-FFLO superconductors. 
Namely, the FFL02 phase occurs only in the presence of strong magnetic field and 
at low temperatures. The phase transition from n to FFL02 phase is discontinuous, 
whereas the transition from FFL02 to the normal state is continuous. 

3. Inter— Site pairing 

In this section we extend our previous study by allowing for the anisotropic <i-wave 
superconductivity. For this sake, we add the nearest-neighbor attraction term to the 
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Hamiltonian (pQ) 

H^H + V £ clciJ^c^ (8) 

where we assume V < 0. For simplicity, we restrict ourselves to FF superconductivity 
on a square lattice. Then, the mean-field Hamiltonian in the momentum space reads 

Hmf = X! ^fco-Cfco-Cfco- 

+ E{[^ff(Q)AS + yrf(fe)A^] c _ fc+QiCfeT 

k 

+ H.c.} - NU cS (Q)\A \ 2 - 2NV\A d \ 2 , (9) 

where A^ = 1/(2 AT) J]fc ^(&)( c -fc+Qi c fcT) andd(fe) = 2(cos cos A; y ). We have assumed 
that the expectation value (c^k+QiCk^) is non-zero only for one particular wave vector 
Q. In clear contrast to case discussed in the preceding section, one has to introduce two 
superconducting order parameters Ao and A^ originating from the on-site pairing and 
the inter-site attraction, respectively. One can straightforwardly calculate the grand 
canonical potential that is of the form 

n = - kT J2 E ln I 1 + exp(-/?£ fejQ )] + £ e_ k+Ql 

a=± k k 

- NU eS (Q)\A \ 2 -2NV\A d \ 2 (10) 
with the quasiparticle energies E kj ± 

Ek,± = of^fcT - £-k+Ql (11) 



2 

± V^feT + ^-fc+Qi) 2 + 4|[/ cff (Q) A + Vd(fc) A 



Numerical results 

The presence of two order parameters Ao and A^ makes the problem numerically much 
more complicated than for the on-site pairing potential only. Though the coexistence of 
the s-wave and c?-wave pairings in the FFLO phase may significantly enhance the upper 
critical field [65], for the sake of simplicity and to avoid too many model parameters we 
restrict our further analysis to the case {7 = 0. 

Generally, the grand canonical potential should be minimized with respect to five 
variables: two components of the wave vector Q, magnitudes of two order parameters, 
and the relative phase <fi between Ao and A^. If the orders do not coexist, i.e., either 
A or Arf vanishes, Q is independent of <fi. However, in order to analyze whether the 
coexistence is possible, we have taken <fi e {0,7r,±7r/2} and minimized Q with respect 
to remaining variational parameters. The resulting phase diagram is presented in Fig. 
|6j Of course, the exact boundaries could be a bit different from those presented in this 
figure if one allows for an arbitrary value <fi. 

Sufficiently strong magnetic field drives the system into the FF state. Depending 
on the values of Q, one can distinguish between two phases marked in Fig. [6] as FFLOl 
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Figure 6. Phase diagram for a square lattice showing the stable superconducting 
phases for various J and h at T = 0. The upper panel shows results obtained for 
V = —2.0t, whereas the lower one shows results for V = —2.51 



and FFL02. In the former case A d ^ 0, Q = (0,Q y ) and Q y <C n. Then, U e s(Q) 
is positive and A = 0. However, in the FFL02 state Q = (ir,Qy), Q y 7r. We 
have found that in the FFL02 phase both the order parameters may simultaneously be 
non-zero. It strongly contrasts with the results obtained in the absence of magnetic 
field, when the system is either in the pure d-wave superconducting state or in the pure 
77-pairing state. Note that for moderate (presumably realistic) values of J and in the 
absence of magnetic the ground state is of purely <i-wave type. 

In order to study the coexisting orders in more detail, we have calculated 
\U e s(Q)Ao\ (see the upper panel in Fig. [7j) and |VArf| (see the lower panel in Fig. 
[7]) in the FFL02 phase. In Fig. [8] we present these data for J = 0.3t together with 
the field dependence of the wave vector Q. Although the dominating contribution to 
the superconducting gap comes from the <i-wave pairing, Ao is non-negligible in the 
FFL02 phase. For sufficiently strong inter-site pairing the boundaries of the FFL02 
phase are almost independent of J, what can be inferred from the lower panel in Fig. [6j 
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Figure 7. |£/ e ff(Q)Ao| (upper panel) and |VAd| (lower panel) for various J and h at 
T — 0. V — —2.5t has been assumed. Dashed lines show show the boundaries of the 
FFLOl and FFL02 phases. In the upper panel \U c s(Q)Ao\ vanishes above the dotted 
line. 




Figure 8. Magnetic field dependence of |J7 cff (Q)A | (a), \VA d \ (b), Q x (c), Q y (d) 
for T = 0, V = -2.5* and J = 0.3*. 



FFLO phase in the presence of pair hopping interaction 



14 



This result can be explained in the following way: increase of the magnetic field shifts 
the wave vector Q and, in this way, modifies U e g(Q). This potential may eventually 
vanish causing U e g(Q)A = (see panel "a" in Fig. ED- Then, the upper critical field is 
determined only by the inter-site pairing. This effect may also be responsible for a non- 
monotonic field dependence of A^. Neither this non-monotonicity nor J-independent 
upper critical field occur for weaker inter-site attraction, shown in the upper panel of 
Fig. [6j Here, both the order parameters are non-zero within the entire FFL02 phase. 

4. Concluding remarks 

Our aim was to investigate the role of the pair hopping interaction for the FFLO 
superconductivity. Probably, this interaction is not a dominating pairing mechanism. 
However, as we have argued in the introduction, in some superconducting systems it 
may become important in the high-field regime whether the FFLO phase is expected. 
Motivated by the pairing symmetry in these systems we have separately studied 
two models where the pair hopping interaction coexists with on-site and inter-site 
attractions. In the former case, the pair hopping interaction lowers the magnetic field 
corresponding to the onset of the FFLO state. In the presence of the inter-site pairing, 
sufficiently strong magnetic field allows for a coexistence of c?-wave and 7/-pairing states 
even though such a coexistence does not occur in the absence of magnetic field. It is 
instructive to compare this result with the recent experimental and theoretical data 
concerning the coexistence of superconductivity and spin-density wave in CeCoIn 5 
[66| 1671 68J. One may formulate a general conjecture, that field-induced breaking of 
the translational invariance of the superconducting phase gives way to other competing 
orders. 

For sufficiently strong pair hopping interaction one may expect 77-pairing state, that 
is robust against the diamagnetic pair breaking. According to our best knowledge this 
phase has not been identified in any known superconducting system. However, upon 
applications of external magnetic field such a system should exhibit FFLO state. In 
contradistinction to the BCS type of pairing, this phase should occur independently of 
the band-width and the orientation of magnetic field. Investigating spatial modulation 
of the superconducting order parameter, one can distinguish, whether the FFLO phase 
originates from BCS or 77-pairing. In the first case the period of modulation is much 
larger than the lattice constant. In the latter case it is of the order of the lattice 
constant and the phase of the order parameter retains its oscillating character typical 
for the //-pairing superconductivity. 
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